Thermodynamic potential with correct asymptotics for PNJL model. 
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I. INTRODUCTION 
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An attempt is made to resolve certain incongruities within the Nambu - Jona-Lasinio (NJL) and 
Polyakov loop extended NJL models (PNJL) which currently are used to extract the thermody- 
namic characteristics of the quark-gluon system. It is argued that the most attractive resolution 
of these incongruities is the possibility to obtain the thermodynamic potential directly from the 
corresponding extremum conditions (gap equations) by integrating them, an integration constant 
being fixed in accordance with the Stefan-Boltzmann law. The advantage of the approach is that 
the regulator is kept finite both in divergent and finite valued integrals at finite temperature and 
chemical potential. The Pauli-Villars regularization is used, although a standard 3D sharp cutoff 
can be applied as well. 

Q \ PACS numbers: ll.10.Wx; 11.30.Rd; 11.30.Qc 

<D 
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One of the most successful and widely used approaches to model the QCD phase diagram in the chiral restoration 
regime is the NJL Lagrangian [l| with the quark condensate as order parameter in the chiral limit, and its extension to 
incorporate the Polyakov loop, the PNJL model where the traced Polyakov loop is an order parameter in the 

O .■ quenched limit and an indicator of deconfinement. For realistic mass values of the light and strange quarks the order 
parameters are approximate quantities, since the chiral symmetry and center symmetry are explicitly broken. The 
J^" 1 . estimates for the transition temperatures from the hadron to the quark gluon phase range from T c ~ 150 — 200 MeV, 
according to recent lattice calculations, see e.g. [22|- (2(|. The nature of the transition as function of the temperature 
and baryonic density, the eventual coincidence of chiral and deconfinement temperatures, the role played by the 
strange quark, the possible existence of a critical endpoint (CEP) and the symmetry breaking pattern associated with 
its location, are some of the topics which are currently under intense study at experimental facilities (GSI, RHIC, 
CERN, LHC), in lattice calculations and effective model approaches, see e.g. [27| - [2!| for recent reviews. 

The NJL model was originally introduced as an interesting field theoretical illustration of the dynamical chiral 
. symmetry breaking phenomena in particle physics. With years it became clear that the effective quark interactions 
capture reasonably well the main features of the low energy meson dynamics. The thermodynamics of the model is 
qq , nowadays extensively used with the hope that some of the basic features of the hot and dense QCD matter are well 
described (at least in the chiral symmetry restoration regime), and this hope has definite grounds (see reviews [30|, 



31]). At the same time, the NJL model thermodynamics, as it now stands, has certain internal problems, related 
with the cutoff dependence of the results. This is a general drawback feature of any non renormalizablc field theory 
model. As a consequence, the asymptotic behavior of some relevant thermodynamic quantities is in conflict with 
some of the well known general laws. The thermodynamic potential separates naturally into the vacuum piece and 
the temperature dependent matter part. It has been reported that the asymptotic description of several quantities 
associated with the thermodynamic potential cannot be achieved using the same regularization criteria on vacuum 
and matter integrals for all observables. For instance, the effective number of quark degrees of freedom has the correct 
asymptotics at large T (Stefan-Boltzmann limit) only when the cutoff is removed from the finite matter contributions 
[32j . On the other hand the removal of the regulator from the matter pieces leads asymptotically to unphysical positive 
condensates of the quarks, except in the chiral limit. In other words the explicit symmetry breaking pattern will be 
restored asymptotically only if the regulator is also kept in the matter integrals. 

In the present contribution we argue that it is possible to achieve consistent regularization criteria for vacuum and 
matter integrals for these observables simultaneously, by obtaining the thermodynamic potential through integration 
of the model gap equations [33| , 34 1 . The integration runs starting from zero over the whole spectrum of the auxiliary 



scalar fields M and is defined up to an integration constant which can be T, fj, dependent. We fix this constant in a 
regulator independent way, by requiring that the low T, [i behavior of the thermodynamic potential yields the correct 
number of degrees of freedom. 

Then we extend this method of obtaining the thermodynamic potential to the PNJL model. With the present 
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treatment we obtain the correct model asymtotics for all observables considered: the number of degrees of freedom, 
the traced Polyakov loop and the dynamical quark masses (or condensates). This is in contrast to the standard 
approach which requires, as in the case of the NJL model, the removal of the regulator from the finite integrals to 
describe the number of degrees of freedom, and its presence for the remaining observables; for example the cutoff 
must be kept to describe the correct asymptotics of the traced Polyakov loop [35|. Although for certain Polyakov 
potentials, of logarithmic form [fij],[3], it is still possible to display its correct asymptotics without a regulator in the 
matter integrals, the physical values of the condensates at high T necessarily require its presence. 

The thermodynamic properties of the NJL and PNJL model have most widely been analyzed regularizing it with 
a 3D sharp cutoff. This has prompted us to use a different regularization to try to clarify the impact that the choice 
of a regulator may have on the asymptotics. We take the Pauli-Villars regularization (PV) [36| which allows to 
check for effects associated with covariance and can be readily used in extensions of the model involving the vector 
mesons, where it is recommendable to use. We take the PV regularization kernel of [37jj , which to our knowledge was 
used for the first time in the NJL model. It has the property that the scalar quadratic and logarithmic divergent 
amplitudes coincide with the ones of the covariant sharp cutoff. Since all the PV regulated integrals can be written 
in 3-momentum space with 3-momentum dependent regulators, it is easy to trace the differences between the two 
methods. This comparison between regularizations will lead us to the conclusion that if one follows the steps indicated 
above to derive the thermodynamic potential, both the PV and 3-D sharp cutoff regularizations can be used to achieve 
the correct asymptotic behavior. 

The paper set up is as follows. After a brief review in section II of the PV regularized fermion one-loop determi- 
nant for the NJL thermodynamic potential obtained in [33| . which we use to introduce the notation,we proceed in 
section III to regularize the PNJL model with the PV regularization. The large T asymptotic behavior of several 
thermodynamic quantities is discussed in section IV. Numerical results for the effective number of degrees of freedom, 
gap equation solutions for the constituent quark masses and Polyakov loop fields are given in Section V, using several 
parametrizations for the Polyakov potential and parameter sets for the many-quark Lagrangian. These results are 
compared with the conventional way of calculating with the 3D cutoff. Conclusions are presented in Section VI. In the 
Appendix C we discuss how the present procedure can be adapted to the 3D sharp cut-off case, leading to compatible 
results for both types of regularization applied on the amplitudes considered. 



II. PV REGULARIZATION OF THE NJL MODEL 



In a recent investigation [33j j we have proven and discussed in great detail how to obtain the correct asymptotics of 
the effective number of quark degrees of freedom at large T for the NJL model, using the PV regularization applied to 
vacuum and matter integrals. This study was done for the NJL model extended to include the 't Hooft 2Nf ( Nf — 
number of flavors) determinant interaction [38| , , [4(| related with the axial U(1)a anomaly, and the eight quark 
interaction Lagrangian 4l| . The latter is proven necessary to render the ground state of the theory stable in the case 
of Nf — 3 [42], 43], 44]. In the text we refer to these extensions of the NJL model as NJLH and NJLH8 respectively 
(whereas the Polyakov extensions are referred to as PNJLH and PNJLH8) 1 . The resulting thermodynamic potential 
as function of the three independent variables Mf = {M u , M d , M s } and T, fx was derived to be [33| 



V(Mf,T,n) = i [-\Ghj 
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In some occasions we will use PNJL to refer collectively to these extensions but we hope it will be clear from the context. 
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to yield the constituent quark masses Mr. Here m/ are the current quark masses, hj = 2 < 0|//|0 >, (/ = 
{u, d, s}), is defined as twice the quark condensate in the symmetry broken phase, and we use the abbreviations 
/i 2 = h\ + h 2 , + h 2 , = h\ + h\ + hj,, and G,K,gi,g2 denote respectively the 49,6(7 and the two 8q coupling 
strengths, of which 171 corresponds to a OZI-violating combination. Given the following specific constraints among the 
couplings, (71 > 0,51 + 3g-2 > 0, G > jiijg) 2 , there exists only a single real solution to the cubic SPA equations [4l| 
and the global stability of the effective potential is guaranteed . For simplicity the case with /i = p u = p<i = p s was 
considered, generalization to account for nonzero isospin chemical potential can be done as e.g. in |45j |. The quark 
one loop integrals in @ are given by: 



r°° 1 
J (M 2 , T, n) = 4 / d\p E \\p E \ 2 p Ai ; E — (1 - n q - n q ) 
Jo tip 



Jo(M 2 ) -4 I d \p E \\p E fp A?E ^pl, (4) 
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f°° 1 /A 2 
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describes the vacuum contribution. The integral J_i(M 2 ,T, p) resulting from the integration over M, 

J_i(M 2 , T, /i) = - / M.M 2 ,T, p)dM 2 = M(M 2 ) + J™ d (M 2 ,T, p), (6) 



splits into a vacuum, 
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and a medium contribution, 
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To obtain the second equality a partial integration has been performed, the boundary term vanishes for the PV 
regularization. The notation n 9 o,n-go, and n q M,n-qM refers to the quark and antiquark occupation numbers for 
massless and massive particles correspondingly, E p (M) — y/ M 2 + p E , E p (0) = \p E \, generically defined as: 

1 _- 1 

1 + e t I + e t 



The Pauli-Villars regulazion kernel that we employ, 



/'V, = Ml-A^)exp(A^), 



corresponds to two subtractions in the integrand |37j . 

The PV result has some unconventional features. First note the occurrence of occupation numbers for massless 
particles in flSJ). These are due to the thermodynamic potential being a direct extension to finite T, p of the effective 
potential which is normalized to V(Mf — 0,0, 0) = and implies the zero mass energy subtraction in the vacuum 
piece (O, (46|, [33j]. As opposed to a 3D sharp cutoff, this subtraction in ([7]), which is a natural consequence of 
requiring the power series expansion of the effective potential to start off at zero when the quark masses vanish, is 
also necessary in the PV regularization scheme to render the 3-momentum integration finite, despite the presence of 
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the regulator. Once one accepts the evidence that thermodynamic consistency cannot be achieved if only the vacuum 
is regularized (shown in the following sections), the massless occupation numbers must be fully considered in the 
regularization of the matter parts. 

Second there occurs the term C(T, /x), with the following significance. The degrees of freedom associated with the 
massless occupation numbers do not vanish in the T — > limit. We use a freedom associated with an integration 
constant associated to the integration of the mass gap equations to impose the boundary condition that as T — > the 
correct limit is obtained. This leads to the M independent quantity 

CP-.,, _ _£ jTdlfcllfcrl^E -> C(T,0) . -?^VT«. (U, 

Since at low T the matter integrals depending on non-vanishing quark masses (and thus also on the PV regulating 
masses) vanish exponentially |33j , C (T, fi) is regularization independent and coincides with the corresponding regulated 
expression (p\p B — > 1). The matter quark-loop part of the thermodynamic potential at finite T, /x, i.e. the </™f d -part 
of the thermodynamic potential, together with the boundary condition term C(T, fi), coincides in the low T regime 
with the matter part obtained with the standard approach when the 3D sharp cutoff is removed (A3 — > 00). 

It is remarkable that the large T asymptotics is as well correctly described. The reason is that at high T the 
regularization of the matter parts becomes relevant. This is taken into account in the PV regulated thermodynamic 
potential (fTJ). In other words we have shown that combining the C(T,fi) term with the regularized expression for 
massless contributions eliminates completely any contribution of massless states to the number of degrees of freedom. 
What remains are the unphysical states associated with the corresponding PV regulating masses which in turn are 
cancelled by the PV regulators of the massive states at leading order of the high T expansion. Thus all auxiliary 
fields introduced in the regulating process are absent in the asymptotic description, leading to the model parameter 
independent and correct result for the number of degrees of freedom. This was proven in all detail in [33| . 

In the present work we adopt a different but equivalent way of looking at these cancellations. We separate the 
term C(T, /x), which does not depend on the regulator, from the remaining regulator dependent matter pieces. Doing 
so we show in section 4 and in the Appendix A that the regularized contributions lead in the high T limit to the 
vacuum expression, up to a sign. The full cancellation of the regularized matter and vacuum pieces at high T allows 
to understand easily the asymptotic behavior of the thermodynamic potential and gap equations. Superficially seen 
this way to group terms may create the illusion that the C (T, /x) is an unphysical Stefan-Boltzman contribution 
associated with massless terms. It is not, since the regularized expression itself contains massless terms which cancel 
it. This is very easily understood by considering the expression (|24[) of the Appendix. The leading terms in the high 
T expansion for the massless and massive terms are equal to each other and to the C(T, /x) term. Their roles can 
be trivially interchanged. This result is the simple consequence that at high T terms of the order M 2 /T 2 can be 
neglected. 

We recall that in the case with 3D cutoff the M — subtractions are not performed and the regulator is a step 
function (p\ 3 — 0(A 3 — |p*e|)). Besides this the integrand of our medium contribution J™ d (M 2 ,T, fi) differs by a 
partial integration from the 3D cutoff one, since the corresponding surface term vanishes in the PV case, as opposed 
to the 3D case. 



III. PV REGULARIZATION FOR THE PNJL MODEL 

We are now ready to extend our analysis to include the Polyakov loop dynamics p|-[Il],|f|,[7|, ( m the present work 
we do not consider the dressed Polyakov loop case, see e.g [ll|). This is effected by coupling the quark and 
gluonic fields in the temporal gauge through the covariant derivative —} = + zA^, with = S^gA^^r 
(A Q , a — 1, ... ,8 are the standard Gell-Mann matrices in color space) and by assuming the temporal component in 
euclidean space-time A4 = %A° to be static and diagonal A4 = A3 A3 + A$\g (the Polyakov gauge). In addition, a 
temperature dependent pure gluonic term is introduced through the Polyakov potential U (<^>, <^>, T), where tj> is the 
traced Polyakov loop and <j> its charge conjugate 

6 = — Tr c L, 1>=— Tr c £t L = Ve^ AxilA \ (12) 

with N c representing the number of colors and V denoting path-ordering. In the framework of the model <f> and 4> are 
treated as classical field variables. The form of the potential U (6, </>, Tj can be determined, for instance, from lattice 
calculations in the pure gauge sector. Several parametrizations have been considered and we shall address them below 
(see Table H|. 
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In the Polyakov gauge it has been shown that the phase of the Polyakov loop appears in the quark action in the 
form of an imaginary chemical potential [9j. This observation allows for a simple extension of the thermodynamic 
potential through the following steps: 



(i) Replacement of the occupation numbers in (j5]) by the tilde quantities 



1 c 

n qM (-E pi jU,T,0,0) = — ^2n qM {E p , n + i (A 4 ) u ,T) 



1 + 3 ((j) + 0e-( £ p-^/ T ) e — t + e 



(13) 



The equivalent expression for hqM {E p , /x, T, </>, </>) can be obtained by replacing /i 4-> — /U and (f> <-> 4>. In absence of the 
gauge field, A4 = 0, one obtains <f> = <f> = 1, and one recovers the occupation numbers ([9]). Note that in the presence 
of a non- vanishing gauge field the limiting value unity for </>, (j) is also reached as T — ► 00, but then obviously the limit 
must be taken simultaneously in all terms of n qi fiq. 

(ii) Specification of the values of <j>, <fi at the integration boundary M = 0. On first sight, it might seem natural to 
systematically implement the substitution (|13|) also for the massless occupation numbers. However this would lead to 
a dependence on <p,<p of C(T, p) which must be field independent in order not to alter the gap equations. Therefore, 
recalling that the boundary M = has been introduced at T = p, = 0, in absence of the Polyakov loop, the values of 
(j), (j) are bound to be unity. 

As a result of these conditions we obtain the final expression for the thermodynamic potential including the Polyakov 
loop as 



V{M f ,T^,^4>) = i (iGh) + nh u h d h s + ^ (h}) 2 + 3g 2 h 
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The minimization of the thermodynamic potential with respect to M u , M4, M s 
equations. 



(14) 

and (j) gives rise to the gap 



IV. ASYMPTOTICS 



At this stage we are already able to draw the following conclusion about the PV regularization on the asymptotic 
value of the total number of degrees of freedom v{T). This quantity at /i = is defined as the pressure difference 
from the zero-temperature value, in Stefan-Boltzmann units 7r 2 T 4 /90, 



v{T) 



P(T)- P (0) 
7r 2 T 4 /90 ' 



(15) 



with p(T) = —V(Mj,(j)*,(j)* 7 T,n — 0), where the starred values Mf,<j>*,<l>* denote the gap equation solutions at a 
given T. In (33|, where the extension to include the Polyakov loop was not considered, we explain how the asymptotic 
behavior of this quantity for large T comes completely and correctly determined by C(T, /j,) independently of the 
model parameters. Now, according to our observation (ii) above, the covariant derivative coupling of the Polyakov 
loop to the quarks does not affect C(T,fi). In the extended version of the model, here considered, once again the 
asymptotic values for the quark number degrees of freedom at large T is determined by C(T, fj.) whereas the Polyakov 
potential U (0*,0*,T) will be the only source of degrees of freedom associated with the Polyakov loop at large T; 
(henceforth we drop the * denoting the gap equation solutions in text and figures, in order to simplify the notation). 
The reason for this is as follows: in the T — > 00 limit we obtain that J_i(M 2 , T, /i, 0, <f) = 0, since then the matter 
part 



J m f(M),T^r 00, n,, 



<1\pe\\pe\ 4 PAp e 



1 



EJM) E P (Q) 



(16) 
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is independent of 4>,4> and equal to the negative of the vacuum expression, J_i(M 2 ), eq. ([7]). The same happens in 
the case of the Jo(M 2 , T, fj,, <j>, <p) integrals (j4}, which reduce to the vacuum expression in the high T limit. Thus the 
hf dependent terms in the pressure also vanish in this limit through the gap equations ^j. 

Therefore we conclude two things. Firstly, since the matter parts become up to a sign identical to their vacuum 



pieces as T — > oo, they need to be regularized to guarantee asymptotic consistency, see also [47], [33(. Secondly, the 
large T asymptotics for the <^>, <j> related number of degrees of freedom is fully driven by the potential term U [<f>, tfi, T) . 
This can also be understood reversing the argument, i.e. leaving the medium part unregularized, which consists in 
discarding the PV subtractions. In this case, using the following representation J3.3J (we consider here for simplicity 
the case with <f> = <j) = 1 and [i = since it does not alter the large asymptotic T dependence of the integral) for a 
term in J™*>d ^ gee ^ e Appendix for a more detailed discussion) 

,,~ ,4 n qM ^ [°° (x 2 + 2xf)t 7^ 4 
dMPEl EjM) ~ Jo dX l + e ^ ~ m T + - 

The last term is the leading contribution in an expansion for T/M >> 1. The unregularized J™^" has then the same 
large T dependence ~ T 4 as the Polyakov loop potential U and can mingle with its asymptotics. It is worthwhile 
pointing out that the complete PV regularized expression, including the C(T, fi) term has also a T 4 dependence 
at large T. However in this case the dependence is due alone to the C(T,fj) term which, as mentioned before, is 
independent of <j) = <j> and is thus decoupled from the large T asymptotics of the Polyakov loop degrees of freedom (see 
also the Appendix). A further remark is in order. The gap equations must be solved together with the SPA equations 
©. Since the Jo(MpT,(j,,<f),(f)) integrals vanish as T — > oo, we obtain that Mf — mf = 0. This guarantees that in 
the large T asymptotics the explicit chiral symmetry breaking pattern of the model is recovered. We emphasize once 
again that this result can only be obtained if the regulator is kept in Dirac and Fermi contributions of the one-loop 
quark integrals. It is easy to understand that failing to do so, the matter parts will diverge at large T, driving instead 
the solutions of the gap equations to Mf — > 0. The same argument as above in (TTT]) can be applied to infer that the 
not regularized J™ 6 " 1 diverges as T 2 . Details are discussed in the Appendix, part B. We shall return to these large T 
asymptotic features in our numerical results presented in the figures below. 



V. NUMERICAL RESULTS 



Before discussing the numerical results we specify in Table Q] several types of Polyakov potentials proposed in the 
literature which we shall employ in our analysis. The Polyakov potential and the pure gauge QCD Lagrangian should 
share the Z(3) center symmetry, spontaneously broken above a critical temperature (for T — > oo the minimum should 
correspond to \cf>\ = 1). The polynomial form U 1 is based on a Ginzburg-Landau ansatz using the terms with the 
required symmetry (the cubic <p terms are needed to break the U(l) symmetry down to the Z(3)), the exponential 
term in U 111 was derived in the strong coupling expansion of the lattice QCD action and the logarithmic term (which 
appears in U 11 , U 111 and U IV ) is inspired by the 577(3) Haar measure of group integration. 

In Table HQ the parameters for the NJLH (sets I [11], II (33[) and NJLH8 model Lagrangians (sets III [H, IV (H) 
are collected, obtained by fitting low lying meson characteristics in the vacuum. Sets I and III have been regularized 
with PV, the remaining with a sharp 3D cutoff, we denote by A generically the PV regulating mass and the 3D cutoff 
in the respective sets. The remaining parameters are defined in the caption of the Table. All parameters are left 
unchanged at finite T and /i, except that in some examples discussed below the regulator will be removed from the 
finite medium contributions. 

Our results are presented in the Figs. 1-5. Fig. la displays the number of effective degrees of freedom v calculated 
for sets I and II of Table 2, in presence of the Polyakov loop. As mentioned before, the PV regularization (full thick 
black line) leads to the correct large T asymptotics, v(T = oo) = 47.5 (dot-dashed horizontal line) determined as the 
sum of quark degrees of freedom, (31.5 for 3 flavors) and Polyakov loop degrees of freedom (16 = 2x8, polarizationx 
color) 2 . If one removes the regulator in the finite matter parts, one obtains the dashed thick black line. As opposed 
to this, the 3-D sharp cutoff shows a large difference between the full regularized expression (full thin curve) and the 
one where only the vacuum part is regularized (thin dashed line, which is almost indistinguishable from the thick 
dashed line). It is worth noting that in the case of the 3D regularization, with or without the Polyakov loop, the 
quark degrees of freedom are largely underestimated above the crossover temperature if all integrals are regularized, 
and only if the finite contributions are not, does one obtain the right order of magnitude. The PV regularization 



2 Note that, as pointed out in 0] this is not the case with U 
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Polyakov potential 


Parameters 


T 4 2 6 V ; 4 V ' 

62(r) _ + ai Z£ + O2 fZ£) 2 + a3 f^ 3 

T \ T J \ T J 


ao ai a2 a3 
6.75 -1.95 2.625 7.44 

b 3 64 
0.75 7.5 


B (0, ^) = 1 - 600 + 4 (0 3 4- 3 ) - 3 (7*) 2 


ao »l a 2 a 3 
3.51 -2.47 15.2 
h 
-1.75 


U 111 U^,T) b , _a__ v 


a b 
664 McV 0.03A 3 


" IV (^- T ) _^< T >^_ & 3 ( ^ +7 3 ) + >4 ( ^ )2 
T 4 2 6 4 y 

- KM -bI </,,</>)] 

24t 2 v ' 


a ai a 2 a 3 
6.75 -1.95 2.625 -7.44 

b 3 fc 4 
0.75 7.5 



TABLE I: Some of the proposed forms for the Polyakov potential: U 1 from @], U" from @, W 7 " from [| and U IV from [l]|. 
Fits to the pure gauge lattice-QCD data result in To = 270 MeV, however in the presence of dynamical quarks a rescaling 
is usually applied and a lower value for this parameter (which appears in U I/II/IV ) is used. The parameter, K, in U IV , 
depends on the choice of parametrization for the quark interactions (in [l5[ it is fitted to reproduce the pressure and transition 
temperature of lattice data). In [6( the parameter b was chosen so as to lead to simultaneous confinement and chiral transitions 
around T » 200MeV (b = 0.03(631.4MeV) 3 ). In [3 U 11 is used with a different parametrization. 



Sets 


m u 

(MeV) 


m s 

(MeV) 


A 

(MeV) 


G 

(GeV- 2 ) 


K 

(GeV" 5 ) 


9l 

(GeV 8 ) 


92 

(GeV -8 ) 


I [48] 


5.3 


170 


920 


8.89 


-687 








II [49] 


5.5 


135.7 


631.4 


9.21 


-740.6 








III [33] 


5.9 


186 


851 


10.92 


-1001 


1000 


-47 


IV [15] 


5.5 


183.468 


637.720 


7.165 


-720.245 


2193 


-589 



TABLE II: Parameter sets I and II are for the NJLH model, using PV and 3D regularizations respectively. Sets III and IV are 
for the NJLH model extended to include the stabilizing eight quark interactions, with strengths gi and g2, with PV and 3D 
regularizations respectively. The remaining parameters are the current quark mass values m u ,m 3 , the cutoff A that indicates 
either the PV regulating mass or the 3D cutoff, the 4q interaction strength G, and the 't Hooft interaction strength k 



fully resolves these problems. The desired approximation to the Stefan-Boltzmann limit is obtained by consistently 
keeping the regulator in all contributions and most importantly by the implementation of the C(T, /j) term (|lip in 
the thermodynamic potential. 

In Figs. 1(b) and 1(c) the constituent quark mass profiles are shown as functions of T for sets I and II, according to 
the different regularizations employed, PV and 3D respectively. In [lj] it has been noticed that the 3D regularization 
applied to the vacuum part only, leads to a sharp fall-off below the current quark mass values, Fig. lc, thin lines. 
This feature is also present if the PV regularization is applied to the vacuum part only, Fig. lb, thin lines. The reason 
was discussed above and related with thermodynamic inconsistency if the matter parts are not regularized. The full 
lines in both Figures correspond to regularizing also the finite matter parts. Only then does one recover in the large 
T limit the explicit breaking pattern of chiral symmetry, in both the PV and 3D regularizations. 

In Fig. 2 we address another issue which is connected with the regularization procedure as well. As mentioned 
earlier, the choice of using a 3D cutoff only in the vacuum integrals was necessary to reproduce the correct behavior 
of thermodynamical quantities in the NJL model without the Polyakov loop [32] . When applied to the PNJL model 



jrocedure results in an asymptotic value for (j) larger 

in the case of 



using the polynomial form of the Polyakov loop potential, this 

than unity when T — > oo (as well as negative susceptibilities) |35|. This can be seen in Figure 2(a) 
the PNJLH model, for vanishing chemical potential. The thick and thin dashed lines correspond to the removal of 
the regulator in the finite parts in the PV and 3D cases respectively. This undesired result is presented as one of the 
reasons behind the adoption of "improved" forms of the potential which include a logarithmic term which diverges 
as <j), <j) — > 1 thus resulting in the correct asymptotic behavior. However we stress that if the matter parts are also 
regularized, (the solid thick and thin lines in 2(a)) the correct asymptotics is obtained. In Figs. (2b) and (2c) we 
show <ft, (ft (full and dot-dashed lines respectively) if a finite chemical potential is introduced, and with all integrals 
regularized, for the PV (thick lines) and 3D (thin lines). Although large differences are observed due to the two 
regularizations at finite T, the asymptotic values will tend all to unity. 
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FIG. 1: Temperature dependence ([T] = GeV) at vanishing quark chemical potential in the PNJLH model with Polyakov 
potential U 1 (see Table HI To = 0.2 GeV) of the following quantities: in Fig. |l(a) is displayed the number of effective degrees 
of freedom, v(T); the thick solid line corresponds to set I from Table ITU with PV regularization applied to vacuum and matter 
part integrals, the thick dashed line is obtained if only the vacuum integrals are PV regularized. Similarly thin lines stand for 
the degrees of freedom calculated with 3D regularization using set II from Table UJl with overall regularization (solid line) and 
if only the vacuum part is regularized (dashed line); the horizontal dot-dashed line represents the Stefan-Boltzmann limit. In 
Fig. 1(b) are shown the dynamical quark masses of the quarks ([Mi] = GeV) corresponding to the set I (PV regularization), the 

. line is 
In Fig. 



thick line refers to the combined use of PV in vacuum and matter parts, the thin line is obtained if only the vacuum part is PV 
regularized; the horizontal dot-dashed lines indicate the current quark masses 
with the same notation as in Fig. |l(b)| for set II (3D regularization). 



1 c 



are displayed the same quantities 




FIG. 2: Temperature dependence ([T] = GeV) of the traced Polyakov loop cf> using the Polyakov potential U 1 from Table UJ 
To = 0.2 GeV (thick black lines correspond to set (I) (PV regularization) and thin lines to set (II) (3D cutoff) from Table ITT1 for 
the following cases: 



2(a) 



/'• 



0, the solid lines correspond to the case where regularization is used for all the contributions 



2 c 



(vacuum and medium parts) whereas dashed lines correspond to the regularization of the vacuum part only; in 2(b)| and 
fi 7^ and all curves are calculated regularizing both vacuum and matter integrals, thick refers as before to PV regularization 
and thin to 3D, and the dot-dashed lines describe the T dependence of (f> (which for finite chemical potential no longer coincides 
with </)); the horizontal dot-dashed line represents the desired asymptotic value for the traced Polyakov loop. 



In the Figs. 3(a) and 4(a) one sees that the introduction of a logarithmic term in the Polyakov potential restricts 
the value of the traced Polyakov loop to below unity independently of the choice of regularization and whether matter 
parts are regularized or not, ensuring the correct asymptotics. However this does not change the conclusion that the 
correct asymptotic values of the constituent quark masses can only be obtained if the finite integrals are regularized 
(see Figs. [3(bJ] ^\ [4(b)] and . 

Finally in Fig. 5 wc show that the effect of inclusion of the 8q Lagrangian for weak coupling strengths g\ and gi does 
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Mi 



M s \ 











(c) 



FIG. 3: The temperature dependence ([T] = GeV, in = 0) of the traced Polyakov loop <j) is represented using the Polyakov 
potential U 11 from Table U To = 0.200 GeV, thick solid lines refer to the parameter set (I) from Table [TT| (overall PV 
regularization) , thin solid lines are for set (II) in 3(a) (overall 3D regularization, they are almost coincident; the dashed curves 



correspond to the regularization of the corresponding vacuum parts only, also almost indistinguishable; the horizontal dot- 
dashed line represents the expected asymptotic value for the traced Polyakov loop. In Figs. 3(b) and 3(c) are displayed in 
thick lines the corresponding dynamical masses of the quarks {[Mi] = GeV), the middle panel for PV regularization, the right 
panel for 3D cutoff, in both cases the current quark mass values (dot-dashed horizontal lines) are approached asymptotically; 
the thin lines denote the cases where only the vacuum parts are regularized, they approach instead zero asymptotically. 




FIG. 4: Same as in the previous Figure (same notation for the curves) using the potential U UI (see Table [fl T = 0.27 GeV): 
thick lines using set I from Table ITU and thin lines using set II from Table HI]. 



not alter the main results discussed so far. This result corroborates our findings 50] 33] that as long as spontaneous 
breaking of chiral symmetry is induced by the 4q coupling strength G, the thermal results of NJLH and NJLH8 models 
will not lead to significant differences. The inclusion of the Polyakov physics does not change this expectation. The 
parameter set IV is taken from [l5;j, where the Polyakov loop was considered together with the NJLH8 Lagrangian, 
and the 3D regularization employed. Note that in this case the minimum of the Polyakov potential U IV at large T 
is at 4>(T oo) = 0.878 and not unity. As discussed before the thermodynamical asymptotics for <fi coincides with 
the minimum of the Polyakov potential if the matter parts are regularized. That's what happens in this case as well 
(solid lines), as opposed to the case where the regulator is removed (dashed lines). The parameter set III taken from 
|33j was chosen also in the weak 8q coupling limit for comparison when the PV is used instead. 
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FIG. 5: The same as in the previous two Figures (same notation for the curves as well) but using the potential U IV (see Table 
HI T = 0.19 GeV and with k = 0.06): thick lines using set (III) from Table \U\ and thin lines using set (IV) from Table EH 



VI. CONCLUDING REMARKS 

We conclude highlighting the main results: 

1- We have presented a way to obtain the thermodynamic potential by integrating the model gap equations and 
making use of an integration constant to ensure the correct low T, // behaviour of the quark number degrees of freedom. 
At the same time the Stefan-Boltzmann limit is correctly described. The quark one loop amplitudes obtained in this 
scheme are calculated using the Pauli-Villars regularization acting both on the divergent integrals related to the Dirac 
vacuum and on the finite integrals over Fermi distributions. A simultaneous description of thermodynamic observables 
related with the NJL models and extensions thereof is accomplished. We addressed three observables: the number 
of degrees of freedom, the quark mass solutions (and implicitly of the condensates) and the traced Polyakov loop 
solutions of the respective gap equations. In particular the large T asymptotic behavior of the number of degrees of 
freedom related with multiquark versions of the NJL model, NJLH and NJLH8, as well as their Polyakov extended 
versions PNJLH and PNJLH8, is correctly described. This is one of our main results. As is well known, in the case of 
the conventional sharp 3D regularization the Stefan-Boltzman limit in NJL and PNJL models is obtained only if the 
matter parts are not regulated. This choice leads however to an incorrect asymptotics for the quark mass and traced 
Polyakov loop gap solutions as indicated in the next point. 

2- We obtain that the constituent quark masses and the Polyakov classical fields (j>, </> calculated with the PV as well 
as the 3D sharp cutoff regularizations lead to similar consistent asymptotic results, when the Fermi distributions are 
regularized. The quark masses approach the current quark mass values at large T. The corresponding condensates 
vanish. As opposed to this a removal of the regulator from the matter parts would result in vanishing quark masses 
in the large T limit, which in turn are associated with unphysical positive condensates. Regarding the <p, <fi fields, we 
show that the Polyakov potential U drives the large T asymptotic behavior of the <j>* , tp gap equation solutions when 
the matter integrals are regulated. Consequently the "overshooting" of these variables at large T above the values 
of the minima of U in the examples shown is consequence of inadequately removing the cutoff from the finite matter 
integrals. 

3- The fact that in point 1 only the PV regularization of vacuum and matter integrals leads to a consistent 
description of the asymptotic behavior, while in point 2 also the 3D regularization applied to vacuum and matter 
parts works, has a very simple reason. The constituent quark masses and <j>, <fi can only detect effects related with 
the gap equations, while the number of degrees of freedom are highly sensitive to the T, fj, dependence of the Mf and 
cj>, <fi independent integration constant C (T, /x) (therefore irrelevant for the gap equations) which is associated with the 
proposed procedure to derive the thermodynamic potential through integration of the gap equations. 

4- Finally, it is possible to extend the method presented here to the 3D cutoff procedure (see Appendix C) and in 
this case both the results of PV and 3D regularizations will be compatible with each other with the same degree of 
thermodynamic consistency for the observables discussed. 
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VII. APPENDIX 

A. Analytical expressions for the large T limit of J™f d (M 2 , T, fj,, <f), <fi) integral in specific cases: 
a) fj, = 0, (f> = <t> = 1; 

The pertinent integral for this case is eq. (fTT|) 



,4 n qM _^ [™^(x 2 +2xf)i 



d| ^ 11 ^ 1 EjW) =T I dx " ■ 



With the variable transformation y = .t + 2.t# it can be rewritten as 



I{M 2 ,T) = T 4 dy Z (19) 

Jo (y/y 2 + a 2 )(l + eVy 2+a2 ) 

with a 2 = ^7^■ The large T limit can be studied by Taylor expanding around a 2 = 0, in complete analogy to [5lj |. 
The leading and subleading orders are easy to obtain, using 

/(M 2 ,T)| a2=0 = ^V (20) 



and 



dI(M 2 ,T) [°° 4 d 

U 2 =o = T / d 2/2T 



da 2 



o dy2 (Vy 2 +a 2 )(l + e^ !/2+Q2 ) 

TT 4 Z" 00 7/ 2 7T 2 

/ ^ V ^_l 2=0 = -T 4 ^, (21) 

2 -/o (Vy 2 + a 2 )(l + eVf 2 +« 2 ) « 



where a partial integration was done to obtain the last integral. The remaining terms 

d 2 I{M 2 ,T) _ 3T 4 f°° dyy 2 d 1 



(da 2 ) 2 2 J dy 2 ^ y2 + a 2)( X + eV V+a 2 ) 



V/ <fr — — _- (22) 



o ( v / y 2 + « 2 )(i + e^ 2+a2 ) 

are identical to the non analytic contribution evaluated in [511 ] . We therefore just quote the final result for the high 
T limit 

I(M 2 , T) = —T 4 - —M 2 T 2 - —M 4 ln[M 2 T- 2 } - —M 4 c + 0(M 6 T- 2 ) (23) 
120 8 32 32 

with c = 27 — § — 2lnir. At M — 0, 1(0, T) = -jfgT 4 . Applying the Pauli-Villars operator to these type of integrals 

appearing in©, we obtain that the regularized integral 

J^l d (M],T — y 00, \x = 0,0 = 1,0 = 1) 

16 rw, f 2 T/„, f 2 , .2 , A 2 " 



■ — [I(Mi,T) — I(M 2 + A , T) + A ——I(Mj + A , T)] 
3 J J aA z J 



+ ^[/(o,r)-/(A 2 ,r) + A 2 ^/(A 2 ,r)] 



1 , M? , M 2 +A 2 „ „ 

f/M^) + A 4 ^(^2— ) + A 2 M 2 ) 



l(M 2 MM 2 )+A 4 ln(l + ^L)). (24) 
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reduces to the negative of the vacuum expression in the T — > oo limit, eq. (0. The large T asymptotics of the 
thermodynamic potential is therefore dictated by the M independent C(T, /x) term. 

As mentioned in section II one sees that in the leading order of the high T expansion terms proportional to 1(0, T) 
and I(M,T) are equal to each other and to the C(T, /x) term. The associated PV regulating masses cancel in the 
leading order as well. 

b) /x = 0,0^ 1,0^ 1 

From lfT3|) one obtains after a variable substitution 

- f°° 1 - 

I[0,0] = / d|p*E||p*E| 4 — — — (h qM (E p ,fj, = 0,T,(j),4>) +nq M {E p ,fi = O,T,0,0)) 

JO Hj p {M) 



OC 



— T Ayy 6 ] ; + 0^0 (25) 

Jo l + e 3 V(^ i ) + 3e 2 v^ i ^ i )0 + 3eV(^ 5 )0 

which is part of the J™f d (Mp T, /x, 0, 0) integral. As before, the first two orders in an expansion around a 2 — lead 
to regular integrals. For the present discussion it is important to realize that the leading term, proportional to T 4 is 
a 0,0 dependent function, for example along = one obtains the closed form expression 

/[0,0]| q2=o = T 4 [^ - 2L U (—— 2 ) - 2L U (— — 2 )], (26) 

360 1 - 30 + V (90 2 - 60 - 3) 1 - 30 - yj (90 2 - 60 - 3) 



where Ln(x) designates the polylogarythm function of 4th order. Therefore the asymptotics of the unregularized 
integral I[<fr, 0], contributes to the number of degrees of freedom associated with the Polyakov loop. This is to be 
contrasted to the regularized J™ (M 2 , T, /x, 0, 0), where in analogy to the example a) the T 4 dependence cancels out, 

being restored through the 0, independent C (T, /x) term. 

The leading order in a 2 of J™ 1 ed (M 2 , T, /x, 0, 0) is the only term that contributes in the limit of vanishing current 
quark masses after the chiral transition has taken place, see Fig. 7 below in part B of the Appendix, and related 
discussion. In this case the results for the gap equations for 0, are known analytically in the T —¥ oo limit and 
serve to illustrate some features of the asymptotics of the more general cases considered in section V. Supposing once 
again that the medium part is not regularized, one obtains in the chiral limit and for vanishing chemical potential, 
for T>T C that 



dl[<t>, 0] ■ - 2T A ^^-^-V^ 2 - 6 ^^ ^^i^^+vW 1 ^)^ (27] 
d<t> N /(90 2 "60-3) ' 1 ' 

with Liz(x) the polylogarythm function of 3rd order, and similarly for the derivative with respect to due to the 
symmetry under f-> of the thermodynamic potential in absence of /x. Since the hight T dependence is the same 
as for the derivatives of the Polyakov potential, one obtains that the asymptotic values of are also steered by the 
medium part of the quark loop integral in presence of the Polyakov loop. In particular as — > 1 the contribution to 

the gap equation for the variable stemming from the quark loop integral is jr ^fe^ 1^=^=1 = ~T~' Thus it is 
evident that 0=1 cannot be the asymptotic solution of the gap equation for for polynomial Polyakov potentials 
that display a minimum at = 1 in the large T limit. 

B. The rcgularization dependence of the quark mass gap equation. 

We illustrate in Fig. 6 for the SU(3) flavor limit the behavior of the regularization dependence of the solutions of 
the quark mass gap equation by plotting the quark one-loop integral Jq(M 2 , T, /i) of eq. ([2]) for different temperatures 
(solid lines) and the T independent quantity 

2iv 2 h 2n 2 r m k , 3. 2 , l9l , , . 

~ Njl - ~ AT [ 7T - (G + 16 h + 4^ + 3 3 ^ ] < (28) 

represented by the dashed line, as function of the condensate h. We use the SPA eqs. ([3]) to express M as function 
of h, since there is a one to one correspondence between these quantities for the model parameters that ensure the 
existence of a stable effective potential. The intersection of the dashed line and solid lines indicates the solution in 
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the dynamically broken phase of chiral symmetry at a certain temperature (here we considered vanishing chemical 
potential and set <j> = <f> = 1). Physical solutions have negative or zero- valued condensates h. In absence of an 
intersection only the trivial solution M = survives (not indicated in the figure, since we are plotting only the non- 
trivial solutions, as we take the ratio -jr?. In the chiral limit (left panel) this quantity is bounded, has a continuous 
behavior and in the limit of vanishing condensate we have lim^o = — ^ ^ 0. The inclusion of a finite quark mass 
however induces a pole (as can be seen in the middle and right panels) which for a physically relevant choice of model 
parameters lies in the positive h side (note that this quantity, expressed as a function of h only depends on the choice 
of the current mass and on the coupling strengths) and the ratio goes to zero with vanishing quark condensate. The 
solid lines correspond to the evaluation Jo for different temperatures. The rightmost panel corresponds to the case 
where the medium part is not regularized, whereas in the other two all contributions are PV regularized. In all of 
them the increase in temperature corresponds to an increase in the negative contribution coming from the medium 
part. From top to bottom the lines correspond to: T = 0,0.1,0.2,0.3,0.4 GeV. 

One can see that in the chiral limit the intersection ceases to exist above a certain critical temperature leaving only 
the trivial solution (M = 0). With the inclusion of a finite current mass there is always an intersection because of the 
divergence near the pole. In the fully regularized version the integral tends to zero with increasing temperature and 
as a result the intersection tends to the solution h = (M = m). When the medium contribution is unregularized 
the integral diverges with the temperature and as a result the intersection tends to the pole (M = 0). 

1.5, 1.5,. 1.5,. 



i.o i.o ; i.o 



-0.20 -0.15 -' 




(a) (b) (c) 



FIG. 6: The intersection of the quark loop integral Jo(M 2 ,T, /i = 0) represented at different temperatures by the solid lines 
(top for T = 0, lower ones for increasing T — 0.1,0.2,0.3,0.4 GeV) with the T independent ratio A given in l)28p (dashed 
line) corresponds to a non-trivial solution of the gap equation (the vertical axis is in GeV 2 units). All curves are plotted as 
functions of the condensate h. Left panel: current quark mass m=0 and overall PV regularization. Middle: m/0 and overall 
PV regularization. Right: m ^ removing the PV regulator from the matter parts. All cases are for the SU(3) flavor limit 
with the parameters of set III from Table [TT] except for the current quark masses which were chosen to be m = MeV and 
m — 20 MeV. See discussion in the text. 

In Fig. 7 we show that in presence of the Polyakov loop and for non-vanishing chemical potential the pattern of 
the solutions for the SU(3) chiral limit with overall PV regularization (left panel) is still similar to the case presented 
on the left of Fig. 6. Removing the cutoff from the matter part (right panel) changes the behavior of the solid curves 
only at somewhat higher temperatures. In this case, since the intersection ceases to exist at low temperatures, i.e. a 
first order transition takes place, the effect of removing the regulator is not noticeable in the solutions. Therefore in 
our study of the regularization effects on the gap equation solutions for the Polyakov loop at high temperature in the 
chiral limit we can start in the phase where the chiral transition has taken place for the quark mass. 



C. The 3D sharp-cutoff case extended. 

The pertinent Feynman amplitudes Jq(M 2 ,T, /i), eq. ((31), and J_i(M 2 , T, /i), eqs. (HJ) , ([7]) , (0 are conveniently 
expressed in terms of the PV regularization kernel (|10p. which can be readily exchanged by other regularization 
kernels. We consider now the case of the 3D sharp cutoff regulator p\ 3 = 1 — Q(\pe\— A3) where the first term leaves 
the integrand unmodified in similarity with the PV regulator. The PV subtractions are then to be replaced by the 
second term of pa 3 ■ This leads to the following results for the integrals considered: 

(i) it is easy to verify that the following vacuum integral which is part of J_i without zero mass subtraction is 
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FIG. 7: The quark loop integral Jo(M 2 ,T, fj,, (f>, (j>) of the PNJL in the chiral limit and the ratio A given in Q28p . using the 
same notation as in the previous Figure and with the same parameters. Here we consider however finite value for the chemical 
potential and the Polyakov loop (fi = 0.2, (j> = 0.4, <j> = 0.6). Left: overall PV regularization. Right: the PV regulator is 
removed from the matter parts. The intersection of solid lines with the dashed line yield a solution to the gap equation for the 
quark mass. See discussion in the text for further details. 



finite, as opposed to the PV case, but non-zero at M = 

d\p E \\p E \ 2 pA 3 E p (M) 

= i(A 3 ^(A§ + A/2)(2A§ + M 2 ) + M 4 ln[M/(A 3 + ^ k\ + A/ 2 )]). 

As such a zero mass subtraction is still required for the effective potential to start off at zero for M = 0. 
(ii) Regarding the medium integral ([5} one has, keeping track of the non vanishing surface term in this case, 



(29) 



8TA2 



Hi 



(l + e (-^))(l + e(-^ i ))) 



(1 + e { -~ t ))(1 + eS-~ t ) 



A\pe\\pe\ A Pk 



i / n qM + HgM _ n q0 + TlgQ 

\ E p (M) ~ E P (Q) 



(30) 



_ 



0, but at T — >• oo it behaves (for p = 0) as 
0(1/T 2 ) and a small dependence on the model parameters of the order of M 2 /A 



The surface term vanishes in the limit T 



-Af + a1Va|+m^- 

-> m 2 /A§ (where m is 



the current quark mass, to which the constituent quark mass approaches asymptotically) remains at leading order of 
the expansion, as opposed to the case of PV regularization, where the boundary term is absent. However this will 
not affect the high T asymptotics, which will be lead by the C(T, /i) ~ T 4 term, as in the PV case, see below. The 
remaining integral behaves qualitatively like in the PV case, with the same leading order behavior in the low and high 
T expansions. To see this, we use again representation (|18[) to obtain (at fi = 0) 



d\p E \\pE\ 4 pA 3 



E P {M) 



dx 



(:r; 2 



■2zf)i 



1 



d.r 



(U+- S /M 2 +Al)/T 



jx 2 +2xf) 
1 + e ^ t 



(31) 



The regulator dependent piece, the second integral, differs from the corresponding expression for the PV subtractions, 
but also vanishes exponentially with T as T — > for non zero values of M and A as in the PV case. Therefore, as 
in the PV case the low T behavior stems again from the same unphysical zero mass occupation states (the regulator 
independent part of the second integral of (|3U|) ); the C(T, /j) term, eq. pT|) . is used analogously to correct it. Therefore 
the low and high T behavior will be the same as in the PV case. 

In Fig. 8 are shown the degrees of freedom for the PNJL model calculated with the advocated way of obtaining the 
thermodynamic potential, using the PV and the 3D sharp cutoff regulators. One sees that the low and high T regimes 
coincide. In the crossover region the number of degrees of freedom is higher for the 3D cutoff as for the PV regulator. 
This is in contrast with what one observes comparing the PV regularization with the conventional 3D treatment and 
the regulator removed from the matter parts, where the behavior is opposite. 
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v 




FIG. 8: Temperature dependence of the effective degrees of freedom using when the regulator is kept in all contributions: 
thick line corresponds to PV regularization, thin line to the usual 3D cutoff procedure and dashed line to the above discussed 
extension of 3D regularization. The choice of parameters is the same as described in Fig. 1. 
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